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The 1996 Brune et al. experiment on vacuum Rabi oscillation is analyzed by means of alternative 
models of atom-reservoir interaction. Agreement with experimental Rabi oscillation data can be 
obtained if one defines jump operators in the dressed-state basis, and takes into account thermal 
fluctuations between dressed states belonging to the same manifold. Such low-frequency transitions 
0^ ' could be ignored in a closed cavity, but the cavity employed in the experiment was open, which 

, justifies our assumption. The cavity quality factor corresponding to the data is approximately 

' Q = 3.31 ■ 10^°, whereas Q reported in the experiment was Q — 7 ■ 10^. The rate of decoherence 

arising from opening of the cavity can be of the same order as an analogous correction coming 
^ , from finite time resolution At (formally equivalent to coUisional decoherence). Peres-Horodecki 

• separability criterion shows that the rate at which the atom-field state approaches a separable state 

' is controlled by fluctuations between dressed states from the same manifold, and not by the rate of 

transitions towards the ground state. In consequence, improving the Q factor we do not improve 
00 ■ the coherence properties of the cavity. 



PACS numbers: 42.50.Lc, 42.50.Dv, 32.80.Ee, 32.80.Qk 



I. INTRODUCTION 



^ , 

^ , In 1996 the group from Laboratoire Kastler-Brossel in Paris reported to date most precise measurements of vacuum 
I ^ i ' Rabi oscillations in cavity QED [l[ . The suggestive plots from [l[ have achieved status of classic textbook illustrations 
of agreement between theory and experiment in quantum optics @, The experimental setup is close to ideal, and 
lO ■ as such is a promising candidate for implementation of quantum logic gates. This brings, however, all the issues one 
K*" [ will encounter in a realistic quantum computer, with qubit-reservoir interactions and various kinds of decoherence in 
the first place. 

When one tries to reconstruct the data from theoretical models one encounters surprising and intriguing difficulties. 
, As far as we know, none of the authors who undertook this task managed to find a convincing explanation of all the 
• ' discrepancies. The main goal of our paper is to discuss theoretical consequences of the fact that the cavity used in 
^ I , the experiment was open. We will see that then at least one new important source of nondissipative decoherence, not 
■ taken into account so far, arises. 

] The first attempt of explaining the data on the basis of master equations was the PhD thesis of Y. T. Chough, 
• • ■ written already a year after the experiment but, unfortunately, not widely known 0, [E]- Chough identified two main 
effects that did not find support in his quantum-trajectory simulations: Too large amplitude damping (probability of 
atomic ground state was too quickly decaying towards 1/2) and too small energy damping (the probability did not 
^ quickly enough approach 1). Chough showed that the data are much better approximated by a semiclassical theory. 
5^ , Theproblem with both kinds of damping was discussed from a general model-independent perspective by Bonifacio et 
al. [6|. Here the conclusion was that the amplitude decay could not be caused by energy losses. The authors localized 
one possible non-dissipative source of the damping, namely experimental time uncertainty At following from atomic 
velocity distribution. Indeed, neglecting energy losses but taking into account uncertainties in the measurement of t, 
one observes that Rabi oscillation of the ground-state probability gets damped towards 1/2. 

Yet another intriguing possibility was noticed by ourselves in [7|. Our original intention was to check if vacuum 
Rabi oscillation experiments are capable of distinguishing between alternative approaches to field quantization. The 
first observation we made was again the problem with too strong amplitude damping. Since in this context this was 
precisely the effect we hoped to see (fields quantized in reducible non-Fock representations of canonical commutation 
relations(CCR) [li, [lOl . [Til [l^ should lead in this situation to vacuum collapses and revivals) we made some effort 
to discriminate between different sources of such behavior. It turned out that the result was inconclusive: The effects 
discussed in [el and ^ can mask one another, unless the uncertainties in At are sufficiently small and one can monitor 
the oscillation long enough to see the revival — if it really occurs. However, the lack of energy decay noticed by 
Chough was clearly seen also in the alternative formalism, so in this respect the situation was not better than in the 
earlier approaches. It became clear that a more reliable theory of dissipation is needed before one can say anything 
concrete about the problem as fundamental as the choice of representation of CCR employed in field quantizaton. 
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Relatively recently Scala et al. [I3 questioned the very basics of master equations employed in the quantum 
optics literature, and proposed an apparently cosmetic modification of the standard formalism. We will describe this 
approach in detail in the following sections, so let us sketch here the main result, which turns out to go back to 
the original ideas of Davies on how to describe system-reservoir interactions in Markovian master equations [13, ITU . 
The observation is very simple but essential for general consistency of the dynamics: The "jump operators" should 
describe jumps between eigenstates of the system Hamiltonian (i.e. dressed states), and not eigenstates of the free-field 
subsystem (free-photon number states) as it happens in practically all approaches employed in quantum optics [T^ . 

We will see that an open-cavity generalization of the model from [13.] we propose allows for reconstruction of 
experimental data even in the idealized case of infinite time resolution. At = 0. The generalized construction involves 
jump operators defined in the dressed-state basis but, as opposed to this from [l^, takes into account also the jumps 
(up and down) between dressed states belonging to the same manifold. Such jumps cannot be excluded in open 
cavities since at T = 0.8K average numbers of appropriate long-wave photons are very large. The cavity of Brune et 
al. was open so the effect has to be included. Of course. At was nonzero and collisions with background gases were 
present. As the latter two effects are mathematically equivalent, we decided to supplement our generalized model by 
the general scheme from [g^. In consequence, we seem to have included the three main ingredients of decoherence: 
Collisions and uncertainty of i, losses of energy, and the external long-wave modes. However, even now we obtain 
agreement with the data only under the assumption that the cavity Q factor was much higher than assumed in 
so some of the problems noticed in [3, ^ do not disappear. A tentative interpretation of this discrepancy is that 
the additional Davies operator we have introduced may arise from system reservoir interactions that are sensitive 
to photon number [l^ (cf. the Appendix). One more observation we make is that theoretical treatments given in 
[1, M, 0, 0] did not properly distinguish between the notions of "true" and "effective" times, somewhat imprecisely 
discussed in and this led to incorrect estimates of decay parameters of the cavity in question. 

The two types of transitions we consider (within a single manifold or between different manifolds of dressed states) 
turn out to have different physical consequences. For example, these are the "ignored" fluctuations within a single 
manifold that control separability properties of the atom-photon system, and this is the rate of this transition that 
plays a role of the damping parameter known from Bloch equation treatments of Rabi oscillations. 

Following the terminology from 13,] we will speak of standard quantum optical "phenomcnological" master equa- 
tions, while the one of Scala et al. will be termed "microscopic". We will confront the data with predictions of both 
models at T = OK and T = 0.8K, and then formulate the open-cavity model at T = 0.8K. In comparison with the 
data we assume exact resonance. Not only will this assumption simplify formulas, but it leads to closed-form solutions 
for a cavity with a Gaussian mode profile. What is interesting, the latter problem can be computed exactly in our 
generalized model, but in the phenomcnological formalism we had to resort to numerical solutions. So the approach 
we regard as more physical, is simultaneously mathematically simpler. 



II. PHENOMENOLOGICAL MODEL AT ZERO TEMPERATURE AND EXACT RESONANCE 



One begins with the master equation 

P = -■i[^,p]+l{apa^ -■^[a''a,p]+), (1) 

where [A, _B]_|_ = AB + BA. The von Neumann part of the dynamics is governed by the Jaynes-Cummings Hamiltonian 
H = hn, where 

^ = y (le)(el-l.g)(5l)+c.oa^a + g(ale)(5l + atlg)(el). (2) 

Creation and annihilation operators a, occurring in Q play also a dual role of Davies jump operators (l^ . [l5| in 
the dissipative part of the master equation. The coupling parameter g is, for simplicity, assumed to be real. Huq is 
the energy separation between the atomic upper ]e) and lower \g) states. The decay parameter 7 is responsible for 
energy losses. 

Let us note already at this stage that this particular model of dissipation assumes that the losses occur only for 
states with non-zero photon numbers: The atom-cavity system loses energy through transitions \g, 1) — > \g,0), but 
]e,0) —^ \g,0) is not allowed. The dissipator in equation ([T]) could be supplemented by an additional term describing 
atomic spontaneous emission at its natural rate unmodified by the cavity. Within the cavity QED regime, however, 
this term can be neglected. In this sense, the presence of the atom inside of the cavity is practically ignored in the 
phenomcnological model. 

In the following discussion we will concentrate on dynamics that starts from the pure state ]e,0)(e,0]. The as- 
sumption implies a convenient constraint on the Hilbert space in question. Indeed, at T = OK the system can make 
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transitions only downwards on the energy ladder, so that the only states we should take into account are: |1) = |e, 0), 
|2) = Iff, 1) and |3) — |(7,0). In this basis Eq. ([1]) is equivalent to the following three sets of coupled first-order 
differential equations 

Pii = «gPi2 - igP2i, (3a) 
912 = -7P22 - igpi2 + «gP2i, (3b) 

Pl2 = -2^12 + «gPll - «gP22, (3c) 
P21 = ~^P21 - jgpii + igP22, (3d) 
P33 = 7P22, (3e) 



Pl3 = -«WoPi3 - igP23, (4a) 

7 
2' 



7 

P23 = -ii^0P23 - ^P23 - «gPl3, (4b) 



and 



psi = iwopai + igP32 (5a) 

7 

P32 = «WoP32 - 2^32 + »gP31- (5b) 

The loss rate 7 does not affect the population pn, while it contributes to the decay of P22 and the increase of P33. 
In other words, the population of the state |e,0) is not directly affected by the dissipation, in contrast to the case of 
\g, 1) (dissipation decreases it) and |.g,0) (dissipation increases it). The coherences pi2 and p2i are also damped but 
with the rate 7/2. 

If we assume that the atom- field system is initially prepared in the state |e,0), the initial condition is pii(O) — 1, 
with all the other elements vanishing. With these conditions the last two sets of equations have solutions identically 
equal to zero, 

Pl3{t) = P3l{t) = P23{t) = P32{t) = 0. (6) 

The remaining five equations imply 

^' 72-16g2 2 ^7" - 16g2 V / 

+ 2(7^:ifg^) -''* {^'^'^^ + e-^*V^) , (7a) 
P22(0 ^ ^^e-** + e-^*V^) - :^^e-i*, (7b) 



72 - 165-2 2 (72 - 16ff2)3/2 

2(72-1652)'' V + J 



7^7 



.-i^fe^^V^^^^-e-^^N/"^"^), (7c) 

(72 - 16<?2)3/2 V J 

7 - 16g2 ^^2 _ 26g2 V y 



UtV^^^ + e-^tV^^^)e-^*, (7d) 
7^ — 16g^ V / 

72 - i6g2 ^72 - i6g2 V y 

. (e^tV^^^ + e-5*V^'^)e-**. (7e) 



72 — 16g^ 



4 



We will be particularly interested in the probabilities p^Qit) = puit) , = P22{t), and P^^o{t) = Pssit). Noting 

that the quantity 7^ — IGg^ appearing in all the above equations takes negative values within the strong coupling 
regime (g is much larger than atomic and cavity relaxation times) we get 





8g^ -it 


I ^ 


16g2 — 7^ 


v/l6g2 - 72 








16g2 — 72 


16g2 — 72 




1 e 
16g2 — 72 


-it 1 7' 
(16g2 - • 



-i*sin^Vl6g^-7^+ ^^^^2 ^2 g"^*cos^Vl6g^-7^ (8a) 
z log — 7 z 



*cos-v/l6g2-72, (8b) 



2 



7 



3/2" "2 

2 J- TR„2 



* sin - ^16g2 — 72 



16g2 ^ 7- ^ ^^^^ - '''^ sin iyiei^. (8c) 

Asymptotically, for t —> 00, we find limf^oo PgV^) ~ ^ ^'^^ limj^ooPg'KO = 0. We conclude that the photon must 
finally escape from the cavity. The solid curves at Fig. 1 show the evolution of p^\{t) and p^\{t) for the atom- field 
coupling characteristic of the Brune et al. experiment The damping parameter 7 is at this stage of our analysis 
chosen arbitrarily, just to show the essential qualitative properties of the dynamics. 

III. PHENOMENOLOGICAL MODEL AT T > — FIRST APPROXIMATION 

Let us now proceed to a more complex and at the same time more realistic case of T > 0. The phcnomcnological 
master equation acquires an additional term related to upwards transitions on the energy ladder, 

^ = +7i(apa^ - ^Ma,p]+) +7T(aVa- ^[aa^p]+)• (9) 

Here, 7^ and 7^ denote rates of transitions downwards and upwards, respectively. Given temperature T, these two 
coefficients are related by the formula (a consequence of the KMS condition [l7| ) 

7T = e~^7x. (10) 

Inserting T and wq as reported in the experiment [1] (T = 0.8K and ujq = 27r 51.099GHz) we obtain 

7l = 671, with e w 0.0466327, (11) 

which, in fact, means that Eq. ([9|) contains a single parameter 7^. The most important formal difference with respect 
to T = is that now the number of states that take part in the evolution becomes infinite, no matter what initial 
condition one takes. This is a consequence of two facts: The form of the thermal density matrix aX t — and the 
possibility of upwards transitions implied by 71 7^ 0. However, for T — 0.8K the contribution to the initial condition 
from all non-vacuum states is not greater than 5% of probability, and the upward transition to higher dressed states 
is of the same order. So the technical assumption we make here is that in the first approximation it is legitimate 
to truncate the Hilbert space to its three-dimensional subspace spanned by |1) — |e,0), |2) — 1) and |3) — |5,0), 
but including the upwards transitions within this subspace. The next approximation would be to include the next 
dressed-state subspace, but within this subspace include only the downwards transitions — this correction will be 
discussed in a separate paper. 

Once we make the simplifying assumption, we have to modify the form of the Davies operators to guarantee that 
higher dimensional subspaces of the Hilbert space will not couple to the subspace in question. Let H be the projector 
on the subspace of interest, 

H = |e,0)(e,0| + |<7,l)(g,l| + |g,0)(<?,0|, (12) 

so that the truncated Davies operators are obtained by the substitution 

a ^ HaH= |5,0)(5,1|, (13) 
at ^ Hatn= |.g,l)(g,0|. (14) 



5 



The modified master equation ([9]) now reads 

P = -i[^,p] 

+7i 0) (5, l\p\g, 1) (<?, 0| - 1) (5, 1|, P]- 

+7T 1) (5, 0|p|ff, 0) (<?, 1| - 0) (5, 0|, . (15) 

As we can see, the dissipative term again does not affect the state |e, 0) of the system. Matrix elements of p between 
basis states split into three groups of equations, 

Pii = «gPi2 - «gP2i, (16a) 

P22 = -7iP22 - «gPl2 + «gP21 + 7TP33, (16b) 
Pl2 = ^Pl2 + ^gPll - «gP22, (16c) 

hi = — 2"'^2i - jgPii + igp22, (16d) 
P33 = 7iP22 -7TP33, (16e) 



7T 

pi3 = -it^oPi3 - igP23 - yPis, (17a) 

7 1 ^ 7t 

P23 = -«W0P23 - jgPl3 7, P23, (17b) 



and 



7t 

Pai = «a;oP3i + «gP32 - Y/331, (18a) 

7 1 ^ 7t 

P32 = «WoP32+«gP31 ^ P32, (18b) 

It is obvious that for 7^ = the equations reduce to ([3a|) -(|5b |l . The situation is analogous to T = 0. Assuming 
that the atom-field system starts from p(0) = |e,0)(e,0| we find that ([TT]) and (|18p again have trivial solutions 
Pizit) — P3i{t) = P23{t) = P32{t) = 0, so that it is sufficient to focus attention on (|16p . Unfortunately, although 
the number of equations is the same as for T = 0, we are not able to determine analytically the eigenvectors of the 
matrix of coefficients appearing at the right hand side of (fT6|) . Numerical analysis shows, however, that predictions 
are basically the same as their counterparts for T — OK, as shown in Fig. 1. Just to have a feel of the role of e for 
the evolution we plot in Fig. 2 the appropriate probabilities for several different es. The range of temperatures we 
consider makes this exercise purely formal, but it clearly shows the tendency for lowering the curves as e grows. 



IV. DESCRIPTION IN DRESSED-STATE BASIS — SCALA ET AL. APPROACH 



The main drawback of the phenomenological model lies in the choice of states determining the possible transitions 
during the energy loss. These are not the eigenstates of the system Hamiltonian H, but those of the free-field 
Hamiltonian hua^a (i.e. I5, 1) and \g,0)). However, in all general theories of system-reservoir interactions one 
assumes that the losses are modeled by transitions between the eigenstates of the Hamiltonian of the system [l3, ■ 
In our case H is the Jaynes-Cummings Hamiltonian representing bound-states of the atom-photon (or spin-oscillator) 
composite system. Atoms and photons are here as inseparable as protons and electrons in atoms and — similarly to 
atoms — should be assumed to loose energy via transitions between appropriate energy levels. The phenomenological 
model is implicitly based on a classical intuition that the system necessarily loses energy by means of losses of light, 
since it is light that interacts with the cavity walls, and not the atom itself (for a detailed discussion of assumptions 
behind the phenomenological model cf. [Tsj). 

An implication for our atom-photon system is that the Davies jump operators a and HaH should be replaced by 
appropriate Davies operators defined in the dressed-state basis. This is the starting point of the "microscopic model" 



6 




FIG. 1: Predictions of the phenomenological model for T = OK (solid) and T = 0.8K (dashed) for g = 477rlO^Hz. The 
dissipation rate for T = is 7 = 0.3g, while 71 — 0.3g and 7^ = 0.04667x for T — 0.8K. Plot (A) shows the probability of 
finding the atom in its ground state, Pg'^(i) = Pg'o(i) + P^\{'t), (B) shows the probability for a photon to be outside of the 
cavity, pf})it). 




0.00002 0.00004 0.00006 0.00008 0.0001 



t[s] 



FIG. 2: Predictions of the phenomenological model for e — 0.01 (solid), e = 0.1 (dashed), e = 0.5 (dot-dashed) and e = 1 
(dotted) and g = 477rlO'^Hz, 71 = 0.3g. As mentioned in the text the restriction to the three dimensional subspace is a gross 
abuse since the corresponding temperatures are, respectively, T ki 0.54K, T « I.IK, T ^ 3.5K and T = 00. 



developed by Scala et al. in [iSj] (further generalizations beyond Markovian and rotating-wave approximations were 
given in [l^lHillSl)- At T = OK the master equation proposed in reads 

p = -i[n,p] + j,(^^\no){n+\p\n+){no\ - ^[\n+){n+\,p]+) +j2{l\no){n^\p\n^){no\ - \[\n^){n^\,p]+){i9) 

where \il±) are the dressed atomic states, |r2o) = IffiO), while 71 and 72 denote rates at which the corresponding 
transitions \^±) \^o) take place. Let us note that fl = H/h is the same as in the phenomenological model: We do 
not modify the system, but only the way it interacts with the reservoir. 

When comparing this equation with ([T]) or ([9]) one notices that it contains two jump operators, |rio)(^^+| and 
|rio)(fi-| expressed in terms of the dressed states, and both of them describe transitions towards the system ground 
state \rto) — l5i 0). Within this model, therefore, both bare states \g, 1) and |e, 0) may decay towards \g, 0) contributing 
to the irreversible loss of energy from the system. In fact, since we are now working with the dressed states which are 
combinations of the bare ones, even a single decay channel produces decays of the two non-ground bare states: |e, 0) 
and \g, 1). Recall, that the essence of the phenomenological model, on the other hand, is that the only state that may 
directly decay is |g, 1). 

From a formal point of view, the essence of microscopic treatments of dissipation is in the way the Davies operators 
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are constructed. In [1^ one gets 

A{^N' ,1 ~ ^N,rn) — \^N ,m) {^N,rn\{a + a^) \Q N> ^i) ^i\ 

= ^Sn.n'-i{VnTT + lmVN)\nN.m){^N+i,i\, (20) 

where l,m — ± and N refers to the total number of excitations in the atom-cavity system. One feature, to which we 
shall refer in the following, is that operators derived this way can describe transitions between neighboring manifolds 
of dressed states, but not within a single manifold. So the model excludes transitions between and |ri_). This 
makes sense in closed cavities since the wavelength of an appropriate transition is of the order of kilometers, while 
the cavity sizes are much smaller (a few centimeters). However, the cavity used by Brune et al. was not closed — we 
shall return to this issue later. 

Assume as before that the initial state of the atom-cavity system is |e,0). In the basis |1) — |2) = and 

1 3) = |r2o) the master equation takes the form 

(21a) 
(21b) 
(21c) 

I Pi2, (21d) 

^) P13, (21e) 
^) P23. (21f) 
As we can see the equation is much simpler than the phenomenological one. The solutions are 

Pii{t) = e-*Vii(0), (22a) 
P22{t) = e-*V22(0), (22b) 
P33(0 = (l-e^*)pii(0) + (l-e-^*)p22(0)-Hp33(0), (22c) 

Pi2(t) = e(-^(^^+-^^-)-^)Vi2(0), (22d) 

Pi3(t) - e(-*(''+-^'")-^)Vi3(0), (22e) 

P2z{t) = e(-*(^^— ^^«)-*)V23(0). (22f) 
In exact resonance the initial condition satisfies 

p(0) - |e,0)(e,0| = i|fi+)(fi+| + i|fi_)(f7_|-i|f7+)(f7_|-i|f7_>(f7+|, (23) 

and the solution can be written as 

p[t) = ie-^*|f7+)(f^+| + ie-^*|f^_)(0_| + (l-ie-^*-ie-^*)|f]o)(f^o| 

_le2'ste-^*|l]_)(0+| - ie-2's*e-^*|f^+)(f7_|. (24) 



Pll = 


71 

-yPii, 




P22 = 


72 

-yP22, 




P33 = 


71 , 
yPii + 


72 

yP22, 


Pl2 = 




- fi_) 


Pl3 




-r!o)- 


P23 = 




-r!o)- 



The corresponding probabilities of the bare states are 



P'^jit) = 1- V^*-ie-^*, (25) 



PTiit) = ie-^* + ie-^*-ie-^*cos2gi. (26) 
4 4 2 



The probability p^"^(t) of finding the atom in the lower state 



= P7^'i*)+P7A(t) = 1 - r"** - ^e-** - ie-^*cos2g< (27) 
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FIG. 3: Predictions of the microscopic approach. Both plots show the probabihty of finding the atom in its ground state under 
different conditions: g = 477rlO^Hz; (A) 71 = O.Ig, 72 = 0, and (B) 71 = O.Ig, 72 = 0.05g. 









< 


2w 





FIG. 4: The experimental setup we consider. The source of atoms (S) and the detector (D) are placed, respectively, directly 
before and after the cavity, so that effects of dissipation outside of the cavity are neglected. The diameter of the cavity equals 
d and duration of the evolution is determined by t = d/v, where v is the atomic velocity. 

is symmetric with respect to 71 and 72 , and it is the sum of these coefficients that controls the rate at which oscillations 
are damped. An interesting situation (population trapping) takes place when one sets 71 = or 72 = 0. In this case, 
the probability of finding the atom in the lower state approaches asymptotically 3/4 instead of 1. In contrast, for 
non-zero 71 and 72 the limiting value is 1, as shown in Fig. [31 

In the following Sections we will try to confront the above two models with the experiment reported in [1,], and 
then propose yet another modification of the master equation, appropriate for an open cavity. 



V. THE EXPERIMENT OF BRUNE ET AL. — ROLE OF GAUSSIAN PROFILE OF THE CAVITY 

FIELD 

The above calculations were performed under the assumption that the atom-field coupling g was constant in time. 
Still, in real experiments the coupling is effectively time dependent since the atoms are moving through the cavity. 
A convenient formulation of the problem is to work in the atomic rest frame, where the moving cavity field can be 
modeled by a time-dependent coupling parameter g{t) [2l|. 

The cavity used in this particular experiment is an open Fabry-Perot resonator consisting of two mirrors of diameter 
d = 50mm each, and with mirror spacing L = 27mm measured on the cavity axis. The mode sustained by the cavity 
has a Gaussian profile characterized by the waist w = 5.96mm. The walls of the cavity (mirrors) are cooled to 
approximately 0.8K which corresponds to the average number n — 0.05 of photons of resonant frequency. The two- 
level system that interacts with the cavity field is implemented by circular Rydberg states of rubidium, with principal 
quantum numbers n = 50 and n = 51, playing the roles of lower \g) and upper |e) levels, respectively. The transition 
frequency between these states is ojq « 2tt 51.099 GHz and the inverse radiative lifetime l/r w 0.03 kHz. The cavity 
relaxation time is T^av = 220/is and g ~ 477rl0^ Hz, hence the parameters fulfill conditions of the strong-coupling 
regime. 
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The experiment is schematically shown in Fig. 4. If the time of flight between the source and the detector is 
the atom arrives at the center of the cavity at t' = t/2. The relation between the distance d and the time t is given 
hy d = vt, where v is the velocity of the atom. The run of the experiment that measures the atomic state at time 
t corresponds to the Hamiltonian H(t') whose coupling is given by g(t') = ge~" ) /™ ^ and the evolution runs 

for < t' < t. Let us note that coupling during the interaction is almost all the time smaller than the maximal 
value g = g{t/2). One expects therefore that the Rabi oscillation (at exact resonance) is not performed with the Rabi 
frequency g but with some smaller effective frequency. Brune et al. assume that this frequency should be computed 
according to 

j^t POO 

g,fft - / g{t')dt' « g / = gV^- = g^^t = gteff. (28) 

Jo J-oo V d 

We will later see that the recipe (|28|) works indeed for system-reservoir interactions involving the dressed-state Davies 
operators, but in fact is not exactly true for the phenomenological model. 

Both the book [23| and the review [2^ make it very clear that the data reported in [l| were plotted as functions 
of the effective time ^cff, and not as those of the true evolution time t. The data obtained in the experiment are 
presented in Fig. O Dots with error bars represent experimental points corresponding to the probability of finding 
the atom in the lower state \g). The solid and dashed curves represent fits by means of the function 

pT'^'M = l-i£^H(l + e~^*°''cos(2gV^M), (29) 

where P{n) is the thermal distribution, and for two different damping parameters 7. 

The fitting formula (j29p has caused a lot of confusion in the literature, and there were two main reasons for that. 
One reason was that the first paper [l| did not clearly state that the time axis involved tcs and not t. The second 
reason can be understood after having solved master equations with time dependent Hamiltonians and our g(t'). We 
will see that inclusion of the Gaussian structure of g(t') effectively changes gi, typical of the constant coupling, into 
gjjffi = gteff, but does not affect the terms ^t. It follows that if we plot the experimental data in the true time the 
fitting function should read 

p-P'(t) = l-if;P(n)(l + e-^*cos(2g,«V^<)). (30) 

But if we decide to rescale the time from t to tcs, as it was done in 1,], then the fitting formula becomes 

pT^'M = 1- 2l]^(")(l + e"^*°''cos(2gV^mieff)). (31) 

In other words, the effective damping parameter is 7/ {-s/TTw/d) and not 7. Many authors, including ourselves [l,[l,[l,0]j 
have noticed that the data presented in [l| seem to correspond to a cavity whose lifetime is 40-50/j,s, and not 220/xs 
as claimed in [ij. It looked like the cavity was five times worse than it was assumed. However, when we insert the 
cavity parameters into j / {^/ttw / d) we find that the effective lifetime of the cavity is around 46/xs, in agreement with 
the data. In Fig. [ntV the time axis corresponds to tcff- The dashed line is the fit with jtcff where 1/7 = 220/is; the 
solid curve is damped by j / {y^w / d)tcff , with the same value of 7. As we can see, an explanation of this discrepancy 
may be trivial. 

One should keep in mind that, in spite of this simple explanation of the "discrepancy", the fitting formulas cannot 
be regarded as theoretical curves, and it is easy to understand why. The long-time asymptotics of both (|30p and (pij) 
is 

hm p-P*(i,ff) = hm p-P*(t) = i, (32) 

whereas we know that the two theories we have discussed so far predict 

limpf(i)= limpf^(t)>i, (33) 

the right-hand-side depending on damping parameters and T, but in all cases the limiting values being much greater 
than 1/2. The predictions of the two models can be mimicked by a different kind of fitting fimction, shown in Fig.03. 



10 




FIG. 5: Probability of the atomic ground state as a function of the effective time ioft . Dots with error bars represent experimental 
points. (A) The dashed line represents (|29[) with 7"^ = 220 • 10~^s; the solid line corresponds to the corrected formula (I31|l . 
hence after unifying the notion of time. (B) We have simply multiplied the probabilities known from the ideal Jaynes-Cummings 
model by exponential factors exp(— 7toff) (dashed line) and exp ( — ^^/j tcs) (solid line). In both cases g = 477rlO^Hz and 
r = 0.8K. 



Physically, the two fits describe exponential damping of two different quantities: The atomic inversion in Fig. [SjA 
and the probability of atomic excited state in Fig. [5j3. It is the later that agrees with theories based on Markovian 
master equations, simultaneously disagreeing with the data. Hence the question remains if one can explain the fitting 
formulas of Brune et al. on a basis of master equations. 

The first question we have to discuss is the precise role of the Gaussian profile of the coupling. We will begin with 
the microscopic model where the form of g^fj will be derived in an exact way. Then we turn to the phenomenological 
model where an exact closed formula will not be available, but the corresponding plot will be produced numerically. 



VI. MODELS WITH GAUSSIAN MODE PROFILE 



Let us begin with a time independent coupling constant g occurring in the master equation p{t) = C(g)p{t). The 
density matrix solution can be formally written as 



p{t) = p/(gH*-*°V(<o). 

Now let g(t) be a discrete n-step approximation of the Gaussian, as shown in Fig. 6. Then 



(34) 



(35) 



where At — t/n, g- — g{tj), tj £ [{j — l)At,jAt], j £ {1, ■ • ■ ,n}. In exact resonance the basis of dressed states is 
independent of g^ . The solution (|22ap - (|22f^ thus depends on coupling constants only via exponents of the off-diagonal 
elements. This allows us to compute the explicit formula valid for the Gaussian, 



Pit) = 



( 



V 



1 -^t 

1 _ Ti+22. 

— — e -1 



-I _ 71+72 

1 -22.4 



tg-2igv^f t 



1 -^t 

1 — e 



1 — e 2" 



(36) 



by means of the limit n — > 00. The coupling parameter g gets indeed replaced by g^f^ = gy/nw/d, as suggested in [l[, 
but the exponents involving 71 and 72 remain unchanged. The probability of the ground state is 



1 T1+T2 t _W 

-e 4 'cos2gV^-t. 



(37) 
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FIG. 6: The interaction time is split here into n = 21 equahy-sized time intervals of length At = t/21. g{t) is expressed in 
units of g and units of time are arbitrary. The distance between the two vertical lines in the central part of the plot is 2w. The 
shape of the Gaussian corresponds to the mode profile inside of the cavity discussed in [ij. In the actual numerical simulation 
of the phenomenological model we use n = 20001. 




FIG. 7: Predictions of the microscopic model at T = 0, (A) with and (B) without the Gaussian structure of the coupling. In 
(A) the atom most of the time experiences coupling much weaker from the maximal value. Parameters are g = 477rl0^ Hz, 
71 = 0.12g, and 72 = O.Olg. t is the true time. 



Fig. [21^ shows how the effective Rabi frequency obtained for a Gaussian coupling differs from the one found for a 
constant g = max4{g(t)} (Fig. [33). 

The case of the phenomenological model turns out to be more complicated, not only because at T > we did not 
manage to solve analytically the underlying eigenvalue equation, but also because even at T = the coupling constants 
enter the solutions (|7a|) - ((7e|) in a very complicated way. So even at T = we would have to rely on numerical analysis. 
The case T = 0.8K is shown in Fig. [51 Having no access to the closed-form solution of the master equation we do 
not know how to switch between true and effective times. So we decided to rescale the experimental data to the true 
time. This is why < < < 430/iS in Fig. [8j3, as opposed to the other plots where < tcs < 91/iS. 

Disagreement between theory and experiment is here at least as bad as in the microscopic model. The structure of 
solutions of these two models does not allow for fitting the data by any choice of the parameters. 



VII. OPEN-CAVITY GENERALIZATION OF THE SCALA ET AL. MODEL 



At T = 0.8K the average number of photons with energy fl^ — f2_ = 2g is n = 354666, but the corresponding 
wavelength is 6378.56 m. If the 5 cm cavity of Brune et al. was a closed one, there would be no reason to take 2g 
transitions into account. However, in open cavities 2j| such large-wavelength oscillations cannot be neglected and 
play a role of a constant noise. An instructive exercise is to compute the temperature that would be needed to produce 
n — 354666 if one took, instead of 2g, the resonant frequency ujq: The result is T = 869770K. The KMS condition 
then implies that the systems jumps up and down between the two states with practically identical rates, Ji ~ 7t- 
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FIG. 8: (A) Predictions of the microscopic model with the Gaussian coupling at T = as a function of t^g. The parameters 
are 71 = 0.12g and 72 — O.Olg. (B) Predictions of the phenomenological model with the n-step approximation (n = 20001) of 
the Gaussian coupling at T = 0.8K, w = 5.96mm, g = 477rlO^Hz. At = t/20001 for each time point (for <t < 430^is with 
1/xs step), t is the true time. The effective Rabi frequency is somewhat smaller than g^g implied by the microscopic model. 
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FIG. 9: Energy levels and decay coefficients used in the generalized model; E± = hQ.±, Eo = hilo. 



The generalized model takes this transition into account but otherwise is the same as the one of Scala et al. The 
new master equation possesses two additional terms responsible for the long-wave open-cavity noise, 

P = -i[^,p] 

+7i(i|f}o>(f^+|p|n+>(r!o| - \[\n+)(n+ipi 
+^a{l\n+){no\p\no){n+\ - ^[\no){no\,p]+ 
+j,[^\no){n^\p\n-){no\ - \[\n-){n-\,p]+) 
+7b(i|f}_)(r!o|p|r!o)(r!_| - i[|^7o)(^^ol,p]+) 
+j3[^\n^){n+\p\Q+){n^\ - l[\Q+){n+ipu) 

+7,(i|r!+)(f}_|p|r!_)(r!+| - (38) 

involving an additional Davies operator |ri_)(r2_|_| and its Hermitian conjugate p^ . The decay parameters are 
explained in Fig. [9l Putting 73 = 7c = we reconstruct the model of Scala et al. To solve psp we employ the 
damping basis method ^25|]. We begin with the time-independent operator £(g)p(t) defined by the right side of 
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and solve its eigenvalue problem. The off-diagonal eigenvectors 



£(g)(|n+)(n_| 

c{g){\n+}{no\ 



71 + 72 + 73 + 7c 



f fn n \ 71 +73 + 7a +76 



|o+)(no|, 



£(g)(|0_)(fio|) = (-i(0_-no) 



72 + 7a + 76 + 7c 



)|0_)(0o|, 



can be supplemented by their Hermitian conjugates corresponding to complex- conjugated eigenvalues. 
To find the next three operators we solve the eigenequation 

£(g)(ar|0+)(^+|+^/|^_)(0_| + ^|^o)(^^o|) = A(x\n+){Q+\+y\Q_){il_\ + z\flo){^o\), 
which is equivalent to 



The solutions are 



Xi 

yi 

Zl 

Ai 

X2 



2 {yic - x{-yi + 73) + 2;7a) 

^(2:73 - y(72 + 7c) + 276) 
i(a;7i+t/72-2(7a + 76)) 



767c + 7a(72+7c), 

737a + 76(71 +73), 
7273 +71(72 + 7c), 

0, 

(71 - 73)(7c - 7a) - (71 + 73)(7l 



Aa;, 
Ay, 
Az. 



72 + 73 - 76 + S), 

2/2 = (71 +73)(73 - 7fe) +73(72 - 7a +7c + S') - 7i76, 

Z2 = -27273 + 71 (71 - 72 + 73 + 7a + 76 - 7c + -S), 

A2 = -;^(7i +72+73+7a+76+7c + 5'), 

X3 = (71 - 73)(7c - 7a) - (71 + 73)(7i - 72 + 73 - 76 - S), 

Vs = (7i+73)(73 - 76)+73(72 - 7a+7c-5')-7i76, 

Z3 = -27273 + 71 (71 - 72 + 73 + 7a + 76 - 7c - '5) 

A3 = -^(71 +72 +73 + 7a +76 + 7c - S'), 



where we have defined 



S 



\J (71 - 72 + 73 - 7n - 76 + 7c)^ + 4(71 - 72) (7a - 7c) 



(39a) 
(39b) 
(39c) 

(40) 

(41a) 

(41b) 
(41c) 



(42a) 
(42b) 

(42c) 
(42d) 
(42e) 
(42f) 
(42g) 

(42h) 

(42i) 
(42j) 
(42k) 

(421) 



(43) 



For the purpose of further reference we list below the set of nine eigenoperators of £(g), defined by the above equations 

(44a) 
(44b) 
(44c) 
(44d) 
(44e) 
(44f) 
(44g) 
(44h) 
(44i) 



Pi 


= xi\n+){n+\- 




.){n.\ + zi\no){no 


P2 


= x2\n+){n+\^ 


^y2\^- 


.){n_\ + z2\no){no 


P-i 


= x3\n+){n+\^ 


^ys]^- 


.){n_\ + z3\no){no 


Pi 


= \^+){^-\, 






P5 


= \^+){^o\, 






P6 








P7 








PS 


= \^o){n+\, 






P9 


= |0o)(n-|. 
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FIG. 10: The dependence of the factor e multiplying 73 in equation (|45cp on temperature. For a wide range of temperatures 
we can simply set 7c — 73; g — 477rlO"^Hz. 



We assume, as before, that the system is in exact resonance. The assumption has one extremely important technical 
implication: The eigenvectors pi, pg, are independent of g. The coupling parameter g enters only into the 
"off-diagonal" eigenvalues A4, . . . , Ag. This property will allow us to derive an exact formula for p{t), even for the 
Gaussian mode profile. 

Up to this point all calculations have been carried out without any assumption about the values of 7a, "ft and 7c. 
This strategy, however, will no longer pay as the resulting equations would be obscure and difficult to analyze. Since 
at r = 0.8K 

7„ = e-^^T^7i w 0.046632771, (45a) 
7ft = 72 « 0.046632872, (45b) 

7c = e~^73 ~ 0.99999773, (45c) 

it is justified to put 7^ = £71, ji, = ej2 with e = 0.0466, and 7c = 73. Fig. [TOl shows that 7c ~ 73 is valid for a wide 
range of temperatures. For temperatures of the order of 10~^K the departures from exact equality amount to O.OO373. 



Rewriting the initial condition p{0) = |e, 0) (e, 0| by means of the eigenvectors we use 

|ri+)(rj+| = — ^ p, ^ 7i+72 + 273 + e(7i+72)-5' ^^ 

2e + 1 7273 + 71 (72 + 73) 2e + 1 4S'(7273 + 71(72 + 73)) 

1 71 + 72 + 273 + e(7i + 72) + -S" 

-P3, (46) 



2e+l 45(7273+71(72+73)) 



2e + 1 7273 + 71(72 + 73) 

7?(e + 1)^ - 5(e7i + 71 + 73) - 73(3e72 + 72 ~ 273) + 71 (-^73 + 73 +72 (e^ - 2e - l)) 

45(7273 + 71(72 + 73))(2e + l)(7ie - 73) 

, -7?(e + 1)^ - 5'(7i + 73 + £71) + 73(72 - 273 + 3e72) + 71 (e73 - 73 + 72 (l + 2e - e^)) 

H ^ ^ Ps- (47) 

45(7273 + 71(72 + 73))(2e + l)(7ie - 73) 
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Accordingly, 



2e + 1 7273 + 71 (72 + 73) 

, 1 1 72(1 + e) - 7172(1 + 3e) + 27173(1 - e) + 273(273 - 726) - S'(7i + 273) 

P2 



-P3 



2 2e+l 45(7273 + 71(72 + 73))(e7i - 73) 

1 1 71^(1 + 7172(1 + 3£) + 27173(1-6) + 273(273 -726) +5(71 + 273) 

2 2e+l 45(7273 + 71(72 + 73))(e7i - 73) 

~Ip4~Ip7, (48) 
^A,p„ X= {1,2,3,4,7}, (49) 



iei 

the latter equality defining the coefficients Ai. Now 

3 



p{t) = ^Aje^^*pj+^4e^^(s)tp4 + ^^e'^^(s)V7. (50) 
j=i 



For t — > 00 the system arrives at a finite-dimensional thermal equilibrium 



limp(t) = -^\n+){n+\ + -—\n^){n^\ + -—\no){no\. (5i) 

t— 2e + 1 Ze + 1 2e + 1 



The coefficients multiplying |r2±)(r2±| are identical as a consequence of our simplifying assumption e 

e ' kT ^' g. The fact that we defined Davies operators in the dressed-state basis results in p(oo) that commutes with 
the system Hamiltonian, as it should. Although wc do not have an exact formula for T > in the phenomcnological 
model, it seems that the asymptotic state will there commute with the free-field Hamiltonian, and thus not with H. 

The asymptotic equilibrium is no longer a Gibbs state, a fact that follows from our truncated-space model. Finite 
dimensional thermal equilibria are known to occur in q-deformed statistics '26] , so probably what we obtain is a state 
belonging to a g-exponential family p7l . [28j , but the problem requires more detailed studies so we leave it open. 

The explicit form of the probability we are interested in is finally 



1 



1 + 2e 



273 — e(7l + 72) — S 71+T2 + 273 + «(ti+72)+5 ^ 273 — e(7l + 72) + S yi+-l2+'^13+'(-ll+12)-S ^ 

^*cos2gi, (52) 

with the asymptotic value 

\im pg{t) = i±l« 0.957. (53) 
t^oo 1 + 2e 

A surprise is that if we put T = 0, 71 = 72 = 0, we find, for the first time in consequence of master-equation theoretical 
analysis, the T — fitting formula of Brune et at. 

P9(i) = ^-^e-^*cos2gi, (54) 

but with the damping parameter whose interpretation is completely different from what was expected. This simple 
test suggests that the open-cavity model may be indeed the correct way towards explanation of the discussed experi- 
ment. Let us remark that formulas analogous to (j54p do occur in (semiclassical) Bloch-equation approaches to Rabi 
oscillations (cf. [11], especially Eq. (11.4) and Fig. 11.2). 

What yet has to be understood is the meaning of the damping parameters — the next Section will be devoted to 
this problem. 

Before we do that, let us briefiy discuss the modification we get for the Gaussian profile. Repeating the argument 
we gave for the Scala et al. model, we find 
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FIG. 11; (A) Dependence ofpg{t) on time and 73. The remaining parameters are g — 477rlO^Hz and 71 = 72 = O.Olg. (B) Pg{t) 
for a longer time (500 /xs) for 73 = 0.2g and other parameters as in the left picture. 73 controls the rate at which oscillations 
decay towards 1/2, whereas 71 and 72 are responsible for the rate at which Pg{t) tends towards 1, reaching it exactly only for 
e = 0. 



pit) = Aie^i* pi + v42e^^* P2 + Ase^'' ps 



and 

— e(7l + 72) — S 71+72 + 2t3 + ^(^1+T2) + S ^ 273 — e(7l + 72) + Jl+-l2 + '^-l3+'(-ll+-l2)-S ^ 

^ 45(2e + l) ' 45(2e+l) 

--e J *cos2gV7r-t. (56) 

The result again agrees with the Brune et al. recipe, with all the restrictions we have made on the meaning of real 
and effective times and parameters. In practice, one also has to include the fact that real measurements involve some 
uncertainty of t in p(t), for example due to collisions of atoms with a background gas and the resulting velocity error. 
An effective density matrix 



PAtit) = / dt'pAt{t,t')p{t'), (57) 

Jo 

where @ 

= At W/At) ' 

satisfies a coUisional nondissipative master equation [2^ [s^l with a double-commutator term —^[Q, [Q, p/\t{t)]]- 
The corresponding ground-state probability 

/"OO 

PgAtit) = / dt' PAtit, t')pg{t') (59) 
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reconstructs the idealized case by limAt^oPg,At(i) =Pg{t)- Inserting ([56)1 into (I59p we find 

1 + e 



l + 2e 

273 - e(7i +12) - S ^ 71 + 72 + 273 + e(7i + 72) + 5 
1 



45(26+1) V ' 4 
273 - e(7i + 72) + 5* A 71 + 72 + 273 + e(7i + 72) - ^ 



4S'(2e+l) 



At 



At 



d2 . 



cos 



At 



arctan 



At 



71+72+273 A^^ , 



(60) 



VIII. CAVITY Q FACTOR VERSUS DAMPING PARAMETERS 



For a constant g the average energy confined in the cavity would be E{t) — hTTflp{t) = hQ{t), where 
070 2e — 1 



fl{t) 



Asymptotically, 



2 2e+ 1 

7ie(a;o + 2g) + 72e(t^o - 2g) + g(7i — 72) — 2a;o73 + S'^o -,'i+-i2+2-ia+'(-ii+~f2)+s ^ 
^ 2S'(2e+l) ^ 

7ie(wo + 2g) + 72e(wo - 2g) + g(7i - 72) - 2wo73 - S'wo -ii+i2+2-,3+'(-ti+~i2}-s ^ 
2S{2^) ' ' ■ ^^'^ 



lhn(f^(t) + ^)=.o(l-^). (62) 



It should be kept in mind that e = does not recover the T = case, since (|6ip has been derived under the simplifying 
assumption that 73 = 7c, so that despite transitions |f2o) l^±) are not allowed, the atom-field system can still 
circulate between and |r2_) with the rate determined by 73. 

An interesting case arises if we assume that 71 = 72 = 7, when r2(t) becomes independent of 73, and reduces to 

t(2£+1) . 



2e + e 



'-t 



m + Y = 26+1 — ^^^^ 

This result can be interpreted as follows. Since 71 = 72, both |ri+) and |ri-) are indistinguishable from the point 
of view of the dissipative processes. In other words, transitions between |ri_) and do not influence the leaks of 
energy from the system. Moreover, at T = the decay simplifies even further 

m + '^ ^ e-i'ioo (64) 

Coming back to (pT|) . with no assumptions about 71, 72, 73 and e, we can compare it to a simple exponential energy 
decay, 

n(t)~n{^ (ji{o)~n{^ (65) 

where Q is the cavity quality factor. 

To find the optimal set of parameters 71 , 72 and 73 we use the iterative Levenberg-Marquardt algorithm for nonlinear 
optimization, setting e = 0.0466 (T = 0.8K) and Q = 7 ■ lO*" as reported in 1]. Straightforward numerical analysis 
shows that the best fit to the right side of (p5|) is obtained for 71 = 72 = 1772IIz and arbitrary 73, which is shown 
in Fig. [T3K. But when inserting these parameters into Pg{t) we obtain poor agreement with the data (Fig. [T3j3; the 
best agreement was found for 73 = 0.057g). If, on the other hand, we try to fit the Rabi oscillation data (Fig. [T5D'). 
the optimal parameters 71 = 72 = 17.73Hz, 73 0.07g imply Q = 3.31 • 10^° (Fig. \WP). The plots in Fig. [13^ and 
Fig. [T2P will not change if we replace in (fCTj) g — 477rlO"^IIz by g = 0. The fits are thus insensitive to the exact 
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FIG. 12: (A) The probability pg, a* (t) given by (|60|l with 73 = and 71 = 72 = 1772Hz, confronted with the Rabi oscillation data 
of [l[. Parameters 71 and 72 are found in the procedure of fitting with exponential energy decay corresponding to Q = 7 • lO'^, 
as reported in [J]. (B) An analogous plots but with Q — 3.31 ■ 10^°, 71 = 72 = 17.73Hz, 73 = 0. t is the true time and 
At = 2.37^s, as estimated in jQ]. 



form of g(t) within the range of parameters typical of the experiment [l|. Moreover, there is no visible change in the 
energy decay plots if we replace ^{t) by 



^Atit) = / dt'pAt{t,t'Mt'), (66) 

with even as large as S/xs, which is greater than the time sampling step employed in the experiment. In consequence, 
the discrepancy between Q we find and the parameter reported in the experiment cannot be explained by finite At. 

Let us now compare the data from [Ij] with ()60l) under the assumption that 73 = 0, that is if the long-wave transition 
|ri+) <-> |r2_) is absent. The estimates made in [6|| led, for the experiment of Brune et ai, to the error At « O.bfis of 
the effective time. For the true time the corresponding value is At « 2.37/xs. Inserting the latter value into (pD|) and 
assuming 73 = we obtain damped Rabi oscillations shown in Fig. [121 The agreement is again good. 

The conclusion is that the effects described in [6] may be of the same order as those arising from opening the cavity. 
Since the cavity was in fact open, it seems that the realistic value of At was probably smaller than the one discussed 




IX. DYNAMICS OF ENTANGLEMENT BETWEEN ATOMIC AND PHOTON DEGREES OF 

FREEDOM 



The last issue we want to address is for how long the atoms and photons in the cavity can be regarded as an 
inseparable bipartite system [sT'] . The problem is important for cavity QED implementations of quantum information 
processing, and is not entirely obvious since we know that there exist two decoherence time scales. Energy losses are 
characterized by 71 and 72 (damping of Pg{t) towards 1), while 73 is responsible for the nondissipative damping of 
Pg{t) towards 1/2. Supplementing the three basis bare states, |10) = |e,0), |01) = \g,l), |00) = |5,0), by the fourth 
state 1 11) = |e, 1) we obtain a bipartite 2x2 composite system whose separability is completely characterized by the 
Peres-Horodecki partial transposition criterion [s^, ■ 

So, let us take the solution p{t), given by (|55p but expressed as a matrix of coefficients evaluated in the basis of 
bare states, and embed it in the space of 4 x 4 matrices, 

/O \ 

Pit) = ° P''''' P''^'' I . (67) 
\ poo,oo/ 
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FIG. 13: (A) The average ^l{t) — r2(oo) of the atom-field system predicted by our model (solid) compared with the exponential 

decay (dots) modeled by (f2(0) — Q,{cx>))e with Q = 7- 10^ as reported in [l|], and 71 = 72 = 1772Hz found by optimization. 
(B) The probability Pg{t) predicted according to 1)56^ with 71 = 72 = 1772Hz, 73 = 0.07g found on the basis of the reported 
Q factor, confronted with the Rabi oscillation data of 0]. Parameters 71 and 72 are found in the procedure of fitting with a 
simple exponential energy decay, while 73 minimizes the sum of squared errors. (C) and (D) show analogous plots, but with 
Q = 3.31 ■ 10^", 71 = 72 = 17.73Hz, 73 — 0.07g; t is the true time. Various combinations of 73 and At lead to practically 
identical plots for Pg,At{t) (e.g. At — 0.5/is, 73 — 0.057g). 
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The Peres-Horodecki criterion states that p{t) is separable if and only if the partially transposed matrix 



/ pio,01 

pio,io 
poi.oi 
V Poi,io 

Poo, 00 



has a non- negative spectrum. It turns out that r{t) corresponding to (|55p has the following eigenvalues 

Ai = {g,l\p{t)\g,l)>0, 
A2 - (e,0|p(0|e,0) >0, 

A3 = i((.9,0|p(t)|g,0) + v/(5,0|p(t)|g,0)2 + 4|(e,0|p(t)|g,l)|2) > 0, 



A4 - - [{g,0\pit)\g,0) - v/(5:0|p(i)|ff,0)2 + 4|(e,0|p(0|ff,l)Pj < 0, 



where 



{e,0\p{t)\g,l) = ie-**(7i- 72 + 273)- 



- 1 



71 - 72 + 73 



2 71+72+^3 j 



sin 2gt. 



(68) 



(69) 
(70) 

(71) 
(72) 



(73) 



The effective measure of entanglement in the system is thus the value of |(e, 0|p(t)|(7, 1)| = |pio,oi(i)l which vanishes 
if and only if 71 — 72 + 273 = and sin 2gt = 0. 
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FIG. 14: Dynamics of separability of p{t). (A) The upper curve represents Pg{t), the same as in Fig. I13D . The lower curve is 
the non-positive eigenvalue A4 of partially transposed p(t). (B) The same curves as in (A) but for a longer time. Entanglement 
decays according to the value of 73 and is unrelated to the time scale of energy losses, controlled by 71 and 72. 

Fig. [HI shows that energy dissipation \^l±) — > |ilo) does not really influence separability properties of atom- photon 
density matrices. The crucial parameter is the one that controls fluctuations between and The fact that it 

is the non-dissipative decoherence that plays the crucial role for quantum computational issues agrees with conclusions 
found by other means in [s^ . One can also think of it as consequence of continuous monitoring of the qubit by the 
noise arriving from outside of the cavity. 



X. CONCLUSIONS 



The formalism where Davies jump operators are defined in the dressed-state basis is mathematically simpler from 
the standard one where the jumps are defined by photonic creation and annihilation operators. Simultaneously, 
when extended to include the long-wave transitions typical of open cavities, the model predicts corrections to Rabi 
oscillations of the same order as those resulting from uncertainties in measurement of t. The so far ignored fluctuations 
\^-) turn out to have fundamental implications for coherence properties of the atom-photon system in cavity 
QED. In particular, these are the transitions that crucially influence separability properties of states in the cavity and, 
hence, the attempts of cavity QED implementations of quantum gates have to minimize both At and 73. Improvements 
of cavity Q factors as well as lowering the cavity temperature yet lower than 0.8K may be less essential. The fact 
that the value of Q we have determined on the basis of the data is some 500 times larger from the value reported in 
the experiment requires further studies. One of the possible interpretations is that the rate at which energy is lost 
in the cavity depends on a^a, and thus for small light intensities the value of Q may be much larger from the value 
one would have found for larger numbers of photons. This interpretation is supported by the observation (l6j that 
the missing Davies operator can be regarded as a consequence of a a^a (g) i? system-reservoir interaction, where R is 
some operator acting in the reservoir Hilbert space. 

We are now in position to apply the above results to the problem of cavity QED tests of different field quantization 
paradigms, in particular to the issue of certain reducible representations of the algebra of quantum fields. There are 
reasons to believe that fields quantized in such representations are much less singular than the standard operator- 
valued distributions, so the question is directly related to various very fundamental issues that have not been subject 
to quantum optical tests as yet. At the level of vacuum Rabi oscillations the reducibility of representations should be 
manifested in terms of collapses and revivals in exact vacuum 7], but observability of the effect — it it really exists 
— will crucially depend on controllability of the coherence losses plaguing cavity QED. These and related issues will 
be discussed in the second part of the present paper [ssj . 
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Appendix: Microscopic derivation of the open-cavity master equation 

The derivation presented in this Appendix was suggested to us by R. Ahcki. Let us start with the Hamiltonian 

H = hfls + ^^B. + hfli, (74) 

where hTls and hTli^ are Hamihonians of the system and the reservoir, respectively. The system-reservoir interaction 
is assumed in the form 

0,1 = A<g)B = |a(a + a1') +/3a'l'a| (K)^gfe(6fc + 4)' C^^) 

k 

(a and (3 are parameters). The term proportional to a is the standard one. The additional part, proportional to /?, 
is the new term proposed by Alicki. Following the procedure described in [l7j we perform the decomposition of the 
operator A: 

A = Y.A{u;), (76) 

UJ 

Ai^) = mAW)- (77) 



e' — e— oj 



n(e) is the projector on the eigensubspace corresponding to the eigenvalue E — he of Hs- Then 

[ns,A{u)\ = -LoA{u), (78) 
[ils,A\oj)\ = +LoA\uj). (79) 
and A^{lo) = A{—uj). Following the standard steps [l3| we arrive at 

p{t) = ^^[Hs,p{t)] + Y.^{u)(A{u:)p{t)A\u)~\[A\oJ)A{u:),p{t)]^ 



w>0 



+ ^7(-^) U{Lu)p{t)A{u^) - ]^[A{u^)A\u^),p(t)\\ , 



(80) 



where we have ignored the energy shifts caused by the environment and assumed 7(0) = (which is true for the usual 
models of thermal reservoirs |17|). 

Now assume that Hs is the Jaynes-Cummings Hamiltonian in exact resonance. The eigenstates of Hs are the 
dressed states 

\^N,±) = -^(|g,7V)±|e,7V-l)). (81) 
The annihilation operators for the atom-field system are given explicitly by 

= (^VN +1+ mm' SN',N+l\i^N,m){^N+l,7n'\ 

+ ^a(^\/N + mm' -Ij 6N'M-l\^N,ni){^N-l,m'\ 

+ ^I3[n + mm' {N - 1)) 6N',N\^N,,n) (f^w,™' | ■ (82) 
To proceed further let us consider the two special cases, = and A^ = 1. For A = 0, 

A{ni,±^no) = -^a\no){ni,±\. (83) 

For A^ = 1 we get 

A{nN',ni' - ^l,ni) = Sn'.2 (y^ + mm' 1^ \ni^jn){n2,m'\ 

+ SN',o\^l.m){^o\ 

+ lf3SN',l\^l,m){^l,m'\- (84) 
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Keeping N = 1 and substituting explicit values for N', m and m' we get the following operators: 



Ain2,+ 


-^1,+) 




+ i)|r!i.H 




A{n2.+ 


- ^1.-) 








Ain2.- 






-l)|f^l.H 


-){^2,- 


A{n2^- 


- ^1.-) 


= 






A{no 






i,±>(f^o|, 




A{ni,+ 












- ^1,-) 


= 






A{ni,^ 




= \mi.A 






A{nu- 


- ^1.-) 




->(^^l,-l• 





(85) 
(86) 
(87) 
(88) 
(89) 

(90) 
(91) 
(92) 
(93) 



In the simplest case T = 0, and with the initial condition we have worked with above, we find 

m - -*[f^s,p(t)] 

+7(f^i,+ - no)a^ Q |f7o) {n+ \p{t)\n+) (f^o I - J i^+lpit)] _ 
+7(f]i,_ - no)a' Q |r!o) {n- \pitm^){no\ - J (f^- 1 , p{t)] _ 
+7(f^i,+ -r!i.-)Y Q|r!_)(r!+|p(<)|r!+)(f}_|-i[|r!+)(fi+|,p(<)]_ 

And this is the equation we have postulated, if 71 = 7(r2i,+ — flo)a^, 72 = 7(r2i,_ — rJo)*^^; a-nd 73 = 
f2i^_)/3^/2. The extension to T > is obvious. 
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